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Past experience with hybrid rockets has shown that certain motor operating conditions are 
conducive to the formation of low frequency pressure oscillations, or flow instabilities, within 
the motor. Both past and present work in the hybrid propulsion community acknowledges 
deficiencies in the understanding of such behavior, though it seems probable that the answer lies 
in an interaction between the flow dynamics and the combustion heat release. Knowledge of the 
fundamental flow dynamics is essential to the basic understanding of the overall stability 
problem. A first step in this direction was a study conducted at NASA Marshall Space Flight 
Center (MSFC), centered around a laboratory-scale two dimensional water flow model of a 
hybrid rocket motor. Principal objectives included: (1) visualization of flow and measurement 
of flow velocity distributions: (2) assessment of the importance of shear layer instabilities in 
driving motor pressure oscillations; (3) determination of the interactions between flow induced 
shear layers with the mainstream flow, the secondary (wall) throughflow, and solid boundaries; 
(4) investigation of the interactions between wall flow oscillations and the mainstream flow 
pressure distribution. The investigation proceeded in two portions: 

The first portion consisted of experimental flow visualization and measurement within a 
test bed developed from a V 3 scale model of an 1 1 inch (28cm) subscale Solid Rocket Combustion 
Simulator (SRCS) hybrid motor which has recently undergone (hot-fire) tests at MSFC. The test 
section was installed in a dual closed-loop water flow facility, originally developed in 1990 to 
investigate fluid flow through porous materials. Oxidizer, or main stream flow, was simulated 
by a water stream entering through injector slots located in the forward section of the "motor"; 
this flow may be injected either axially (along the motor axis) or radially. Fuel grain "burning" 
(mass injection only) was simulated by a secondary water stream entering the test section through 
100p sintered bronze porous plate material. The test section side walls were constructed of clear 
acrylic for flow visualization and optical velocity measurements. Flow visualization was 
accomplished through injection of very small helium bubbles into the mainstream and/or 
secondary flow streams and recording their movement with a 1000 frames/second video recorder. 
Qualitative flow field mapping was accomplished with an LDV system. A test section schematic 
is shown in Figure 1 . A complete description of the system setup and results obtained is given 
in Reference 1. 

The second portion of the investigation consisted of the generation of an analytical model 
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to supplement the experimental data obtained, and is the subject of this report. The model is an 
extension to hybrid rocket motors of an analysis originally developed by Ciucci and Jenkins for 
solid rocket motors 2,3 ' 4 ' 5 . The analysis will handle unsteady, two-dimensional (or axisymmetric), 
incompressible or compressible flows characterized by the presence of one or more jets issuing 
into a finite volume. Boundary walls may be porous or non-porous. The two-dimensional, 
unsteady, compressible, full Navier-Stokes equations are employed. The overall computational 
space is characterized by one or more planes of symmetry, with the analysis restricted to a single 
plane of symmetry for computational efficiency. Though not applicable to the present application 
(incompressible flow), the model can easily be extended to incorporate locally imbedded 
supersonic regions. The principal governing fluid flow equations are as follows: 

Conservation of mass: 
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Conservation of Momentum: 
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Conservation of Energy: 


dpe + dpue + d£ve _ 
dt dx dy 


dx dx dy dy 


du + dv, 
dx dy 




du 

dy 


dv^ 

dx' 




In the above equations t is time, p is density, p is pressure, u and v are velocity components 
associated with x andy spacial directions, respectively, T is temperature, e is energy, and p and 
k are flow viscosity and thermal conductivity, respectively. The <I> terms are source terms whose 
exact form depends upon whether the computational domain is two-dimensional or axisymetric 
(2-D in this case). Though a turbulence model can be built in to the analysis, only laminar 
flows are considered here. 


Numerical Technique 

The numerical solution of the equations of motion is obtained utilizing an explicit, 
time-dependent, predictor-corrector finite difference method developed by MacCormack 6 . In 
order to damp numerical oscillations induced by the gradients in the flow field associated with 
the developing jet (and wall) flows, a fourth-order damping scheme introduced by Holst 7 and 
modified by Berman 8 and Kuruvila 9 is utilized. This scheme involves certain free adjustable 
parameters, C x and C y , usually referred to as damping or dissipation coefficients. Generally, C x , 
C y should be such that 0 <C X ,C y < 0.5. General nomenclature for the computational domain 
is shown in Figure 2. 


Computational Grid 

A uniform rectangular grid is utilized in the upstream plenum section, extending a 
distance of 1.5 times the CP height into the CP portion of the domain. For the extended portion 
of the domain, which encompasses the outflow region, a scheme introduced by Cebeci and Smith 
is employed 10 , whereby grid spacing is increased by a fixed percentage from an initial value. 
General nomenclature for the mesh generated in this manner is shown in Figure 3. 


Initial and Boundary Conditions 

At time t = 0 velocity components u and v are everywhere set to u=v=0, and pressure and 
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HYBRID MOTOR COLDFLOW MODEL 
GENERAL NOMENCLATURE 






Figure 2: General Nomenclature, Computational Domain 
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HYBRID MOTOR COLDFLOW MODEL 
GRID GENERATION NOMENCLATURE 


(NX1P.LY-1) (LX4.LY-1) 



(1,1) (LX.1) 


Figure 3: Nomenclature, Computational Grid 
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temperature are set to ambient values everywhere in the flow field. Density and internal energy 
are obtained from the equation of state. Boundary conditions must also be specified for all of 
the dependent variables, u, v, p, e, along with corresponding values of p and T. Along all solid 
walls of the computational domain a "no slip" condition is enforced, that is, u=v=0 for all time 
t. An adiabatic wall condition is imposed at all solid walls, as is the condition that the normal 
pressure gradient equals zero there. This latter condition is obtained from the momentum 
equation applied in the direction normal to the wall, as in a boundary layer type flow situation. 

The top boundary of the computational domain is an axis of symmetry. Special note 
should be taken of the fact that, unless the outflow boundary is properly chosen, the downstream 
boundary conditions are not known. The approach adopted here is to continuously move the 
downstream boundary as the solution proceeds in time so that the initial disturbances (pressure 
waves) never reach it. Thus, it is possible to assume unperturbated conditions at the outflow 
boundary, requiring that both velocity components are zero, u=v=0 and that pressure, temperature 
and density are equal to their corresponding initial, ambient values. 

The remaining boundary conditions are for the jet (oxidizer) flow and for the porous wall 
flow(s). Both the upstream plenum wall and the CP wall are porous, as in the experimental test 
rig. At time t=0, all jet and porous wall flow rates are zero. These flows are then ramped up 
linearly with time up to specified values, at which point they remain constant. Thus, the flow 
field variables at the inflow boundary do not initially remain constant with respect to time at the 
inflow boundary; this period is, however, very short as compared to the overall calculation 
interval. 

Computer Code and Sample Results 

The computer code (Appendix A) is written in Fortran-77 and is currently run on a Sun 
Systems Sparc 10 workstation. Simulation of the (incompressible) flow in the hybrid model test- 
bed was accomplished through Reynolds number scaling of the measured flowfield velocities to 
a comparable flow in a gas having the properties of air at "room" temperature. The resulting 
velocities were well within the incompressible flow regime (Mach number < 0.2). All input units 
are SI, but provision is made for US units if the user so desires. Required input is documented 
by annotation within the code listing (Appendix A) and will not be discussed further here, except 
that geometrical input nomenclature for a typical (spacerless) configuration is shown in Figure 
4; the corresponding calculation stations are shown in Figure 5. A sample listing of typical input 
for the case of a radial injector configuration with a spacer is shown in Appendix B. 

Figure 6 illustrates a comparison of experimental results" with computed results for axial 
flow injection through a non-porous wall. Figure 7 provides a similar comparison (same 
reference) for a radial flow injection pattern. Figures 8 and 9 compare LDV data (from portion 
1 of this investigation) and calculated data, respectively, for the case of axial flow injection with 
no spacer (inhibitor, in the actual motor). Figures 10 and 11 provide the same comparison for 
a radial injection pattern. Figures 12-15 repeat the comparisons for the case where a restrictive 
spacer is present between the upstream plenum and the CP portion of the model. 
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HYBRID MOTOR COLDFLOW MODEL 
INPUT NOMENCLA TURE 


X 1 (NX11) X 2 (NX1) X 4 (NX2) 



MLENGTH = OVERALL MOTOR LENGTH 

NYIG = NO. GRID POINTS IN INJECTOR EXIT PLANE 


Figure 4: Computer Code (Geometrical) Input Nomenclature 
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HYBRID MOTOR COLDFLOW MODEL 
CALCULATION NOMENCLATURE 


(NX1P.2) 


(LX4.2) 


(NX1P.NY2PP) 
(1.NY2PP) 



(LX4.LMY) 


(1,LY) 


(LX.LY) 


Figure 5: Computer Code (Geometrical) Calculation Nomenclature 
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Figure 6: Comparison of Code Results with Experimental Data 11 
Axial Flow Injection (No Wall Blowing) 
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Figure 7: Comparison of Code Results with Experimental Data 11 
Radial Flow Injection (No Wall Blowing) 
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DV TEST DATA: AXIAL FLOW INJECTOR 

(No Inhibitor) 



Figure 8: LDV Data, Axial Flow Injection, No Inhibitor 
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AXIAL FLOW INJECTOR: NO INHIBITOR 
(CALCULATED FLOW) 



Figure 9: Calculated Data, Axial Flow Injection, No Inhibitor 
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RADIAL FLOW INJECTOR: NO INHIBITOR 



Figure 10: LDV Data, Radial Flow Injection, No Inhibitor 


14 





Figure 11: Calculated Data, Radial Flow Injection, No Inhibitor 
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AXIAL FLOW INJECTOR: WITH INHIBITOR 



Figure 12: LDV Data, Axial Flow Injection, Inhibitor Present 


16 



IAL FLOW INJECTOR w/INHIBITOR 
(Calculated Flow) 



Figure 13: Calculated Data, Axial Flow Injection, Inhibitor Present 
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Figure 14: LDV Data, Radial Flow Injection, Inhibitor Present 
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ADIAL FLOW INJECTOR w/INHIBITOR 
(Calculated Flow) 



Figure 15: Calculated Data, Radial Flow Injection, Inhibitor Present 
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APPENDIX A 


Computer Code Listing 
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APPENDIX B 


Typical Code Input 



0.2 0.35 0.35 0 18000 
.2700 100 . 20.0 6.0 

.4 287 . 29 . 0.0001 0.007 

01.353447 298 . 

.0127 0.0381 0.020108 0.01905 0.1016 0.0137584 

.047244 12 36 96 19 18 13 

.05 

00 . 1.177 287 . 0.0261 0.0000185 
6 1 0.01016 7810 



